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Abstract 

Anomalous dimensions of high-twist Wilson operators in generic gauge theories occupy a band 
of width growing logarithmically with their conformal spin. We perform a systematic study of 
its fine structure in the autonomous SL(2) subsector of the dilatation operator of planar M = 4 
super Yang-Mills theory which is believed to be integrable to all orders in 't Hooft coupling. We 
resort in our study on the framework of the Baxter equation to unravel the properties of the 
ground state trajectory and the excited trajectories in the spectrum. We use two complimentary 
approaches in our analysis based on the asymptotic solution of the Baxter equation and on the 
semiclassical expansion to work out the leading asymptotic expression for the trajectories in 
the upper and lower part of the band and to find how they are modified by the perturbative 
corrections. 
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1 Introduction 



Wilson operator product expansions (OPE) [T] are the main tools in analyses of light-cone 
dominated processes, with deeply inelastic scattering being the most prominent representative. 
These expansions underpin the incoherence of strong-interaction phenomena occurring at dif- 
ferent space-time scales and yield separation of the long- and short-distance dynamics of the 
underlying theory of strong interactions. The former is encoded in hadronic matrix elements 
of composite Wilson operators built from elementary QCD fields X = F^, if), x/j with properly 
contracted color indices to form SU(3) singlets and an arbitrary number of covariant derivatives 
D^. Among these operators there exists a subclass of the so-called quasipartonic operators which 
play a distinguished role and which we shall study in this paper. These operators arise in the 
OPE expansion of light-ray operators of the following (schematic) form 

X( Zl n)X(z 2 n) . . . X(z L n) = O^, • • • , ^i)O^(O) + . . . (1.1) 

7V>0 t 

where is the light-like vector, n 2 = 0, while the scalar variables Zj define light-cone coordinates 
of the fields and the coefficient function Cl,n( z u, ■ ■ ■ > ^li) is a homogenous polynomial in = 
Zi — Zj of degree N. It is tacitly assumed that the gauge invariance of the nonlocal operator 
on the left-hand side is restored by inserting Wilson gauge links between the fields and by 
appropriately taking traces over the fields' color indices in corresponding SU (N c ) representations. 
The quasipartonic operators O^ L (0) are built from L quantum fields (not necessary identical) 
and N light-cone derivatives acting on them. They belong to the SL(2) sector of the conformal 
group and the parameters L and N define their twist and Lorentz spin, respectively. Conformal 
symmetry allows one to separate quasipartonic operators into towers with total derivatives like 
(n ■ d) p O^ L (0). The latter operators are shown by ellipses in the right-hand side of ( 11. ip . It is 
also known that starting from twist L = 3, one can construct several operators with the same 
quantum numbers N and L. The sum over I in the right-hand side of ( 11. ip is meant to enumerate 
such operators. The physical meaning of this quantum number will become clear in a moment. 

Due to the interaction between partons, e.g., quarks and gluons, the Wilson operators mix 
under renormalization. Diagonalizing the corresponding mixing matrix we can construct the 
operators O$ L (0) in such a way that they have an autonomous scale dependence governed by a 

set of anomalous dimensions 7^ L . The matrix elements of the Wilson operators O^ L (0) encode 
information about partonic structure of hadrons interacting with hard probes, i.e., photons and 
weak bosons, and their anomalous dimensions control the scale dependence of various physical 
quantities. Leading contributions to physical cross sections of high-energy scattering are associ- 
ated with twist-two operators [L = 2), while power-suppressed effects with higher twists (L > 3). 
In particular, the matrix elements of twist-L operators between the vacuum and a hadron state 
determine partonic distribution amplitudes [2] 

(0|O^(0)|P) M ~ / [dx] L ^ ) (x 1 ,...,x L )f(x 1 ,...,x L]f i), (1.2) 

where Xi are the momentum fractions of partons in the hadron with momentum P^. Here, the 
integration measure includes the momentum conservation constraint [dx\i = dx\ ... (1x^8(1 — 
Ylt=i an d the weight function (p^(xi, . . . , x£) is uniquely defined by the coefficient function 
Cjv^ in Eq. (II. ip . Both sides of this relation depend on the renormalization scale \x and the 
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Callan-Symanzik equation for the Wilson operators can be translated into the evolution equation 
for the distribution amplitude f{x\, . . . ,2^;//). Solving the evolution equation one finds that 
f(xi, . . . , Xl', [J>) can be decomposed into the sum of different components enumerated by I each 
having an autonomous scale dependenc^] governed by the anomalous dimension 7^ £ 

f(x u ...,x L ]fi) = V V0$ L (xi,...,x L )( — ) UnIlM, (1.3) 

where the functions 4>^ L (xi, . . . , Xl) are orthonormal to the set of the weight functions in (jl.2p . 

fio is some reference scale and a$ L (no) = (0|O^ L (0)|P) M0 is a nonperturbative parameter nor- 
malized at /i . 

All ingredients of the distribution amplitude in the right-hand side of (11.31) are important for 
successful phenomenological description of experimental data. Theoretically, however, these are 
the anomalous dimensions which represent "easier" observables for analytical treatment since 
they can be calculated starting from perturbation theory in gauge coupling g YM . The anomalous 
dimensions J^ L are defined as eigenvalues of the mixing matrix which in its turn represents the 
dilatation operator of the underlying gauge theory in the SL(2\ M) sector. For given twist L, 
the mixing matrix can be interpreted as L-particle quantum mechanical Hamiltonian and the 
problem of finding 7^ L can be reduced to solving the eigenproblem for this Hamiltonian. The 
corresponding energy levels are parameterized by the quantum numbers N and I. It is convenient 
to organize the energy spectrum into a collection of trajectories labelled by t. Each trajectory 
defines a smooth functions of N whose value coincides with 7^ for N being nonnegative integer. 
As in quantum mechanics, the energy levels of the Hamiltonian do not cross and, therefore, the 
trajectories with different i and £' do not cross either. This allows us to order them for each L 
and arbitrary iV > as 

1 { Z<1n,l<--^ £ = fixed, N>0, (1.4) 

where 7^ is the ground state trajectory, 7^ is the first excited trajectory and so on. We would 

like to stress that the relation (11 .4]) holds to any loop order. In other words, if 7^ < 7]^ L at one 
loop, then the same should be true for arbitrary coupling. Then, it follows from ( 11.31) . that for 
\i — ► 00 the sum over components in the right-hand side of ( 11.31) is dominated by the contribution 
of the ground trajectory i = while the effect of the excited trajectories on the observable is 
suppressed by [i~ Ae with = 7^ — 7^ > 0. This explains why the ground state trajectory 
plays a special role in the analysis of the distribution amplitude (jl.3p . 

More than this, anomalous dimensions and the properties of the trajectories (11.41) are of 
interest in their own right as they were found to reflect a (hidden) symmetry of the underlying 
gauge theory. In particular, the one-loop QCD dilatation operator, whose eigenvalues determine 
one-loop anomalous dimensions of conformal Wilson operators, was found to possess integrability 
symmetry in the so-called aligned-helicity sector. It can be made manifest through the map of the 
dilatation operator to the Hamiltonian of an exactly solvable SL(2) Heisenberg magnet [31 HI EJ [6] . 
This integrable structure was naturally embedded into a more general SU(2, 2|4) spin chain [3,(8] 
spawned by the M = 4 super Yang-Mills theory which encodes anomalous dimensions of all single 



1 Strictly speaking, this relation only holds in conformal theory. In theory with nonvanishing beta-function, 
the /u-dependence is more complicated. 
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trace operators of this superconformal theory. The one-loop consideration was extended to all 
orders in 't Hooft coupling g 2 = g 2 M N c / (47T 2 ) P, [H)] thus providing a framework for strong 
coupling analysis of anomalous dimensions, on the one hand, and a test of the gauge/string 
duality between maximally supersymmetric SU(N C ) Yang-Mills theory in the large- iV c limit and 
the type IIB string theory on AdSsxS 5 background [11] which relates anomalous dimensions of 
Wilson operators to the energies of corresponding string configurations [T2l fl~3] . 

For the twist-two operators, the eigenvalues of the mixing matrix reside on a single trajectory 
which is a function of the conformal spin N only. At large conformal spins the anomalous 
dimensions displays logarithmic Sudakov scaling [TJl [15] , 

lN,L=2{g) = 2T cnsp (g) In AT + O(N ) , (1.5) 

with the overall coefficient determined by the cusp anomalous dimension. For higher twists L > 2, 
the analytical structure of anomalous dimensions is more complex. In this case the anomalous 
dimensions are not specified by the conformal spin alone and form families of nonintersecting 
trajectories (11.41) described by the continuous functions 7^ L which depend on the 't Hooft cou- 
pling g. At large N the trajectories occupy a band of width which grows logarithmically with N 

mm- 

2T cusp (g) InN < ^ L (g) < LT cnsp (g) InN. (1.6) 

As one approaches the strong coupling regime, one anticipates that these trajectories do not 
intercept however their distribution inside the band gets modified. The minimal trajectory has 
the same leading large N behavior as for the twist two [161 E] 

7 $ L (0) = 2r cusp (s) IniV + O(N ) (1.7) 

but subleading O(N ) terms are different as compared to (11. 5p . 

We would like to stress that the relations (II. 51) and ( 11. 61) are valid in any gauge theory, even 
in the 5*L(2;M) sectors which are nonintegrable. The value of the cusp anomalous dimension 
depends however on the particle content. In generic (supersymmetric) Yang- Mills theory, the 
cusp anomalous dimension is known to two loops [16] , whereas in M = 4 SYM theory it was 
explicitly calculated in the first four orders of perturbation theory [161 O EE 1201 [21] . Later in 
the paper we will use its expression to three loops 

Ws) = s 2fcr W-i = g 2 - fr*g A + i^y + o( g 8 ) . (1.8) 
fc=i 

The integrability allows one to describe the fine structure of the trajectories. To one-loop 
order, the anomalous dimensions in the integrable SL(2) sector coincide with the energy spectrum 
of the Heisenberg SL(2) spin chain of length L and the total spin N and can be found from the 
Bethe Ansatz. The Bethe roots take real values only and condense at large iV on the union of 
intervals on the real axis [<7i , 02] U [as , 04] U • • • U [02L-3 , 022,-2] . The values of <7i can be determined 
by solving the Bethe Ansatz equations in the semiclassical approximation. The corresponding 
solutions are parameterized by the set of integer numbers 

L-l 

£ = {£i,e 2 ,-^jL-i}, £i>o, J2 e * = N > (°) 

i=l 
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where li counts the number of Bethe roots on the interval [021-1, <72i] with i = 1, . . . ,L — 1. The 

(£) 

same set of integer numbers determines the anomalous dimension ^ NL in the semiclassical ap- 
proach. In order to describe a particular trajectory in the spectrum of the anomalous dimensions 
7^ L , it suffices to fix L — 2 independent integers £\, . . . , Il-2 and allow iV to take positive con- 
tinuous values. Going from the lowest to the highest trajectory inside the band (11.61) amounts to 
counting all integers satisfying (11. 9p . The one-loop expressions for the trajectories have been 
found in [22j [4], [U [6] and sum rules for the excited trajectories to high loops were discussed in 
[23]. 

The properties of the trajectories depend on both the spin N and twist L. At large N the 
leading asymptotic behavior of the anomalous dimension jjy L ~ IniV does not depend on the 

twist L. The subleading corrections to ^ N L cease to be twist-independent and, as was shown in 
Ref. [IT], their properties depend on a hidden parameter £ = -^-lniV. Namely, for L, N — ► oo 
and £ 3> 1 the minimal anomalous dimension belonging to the ground state trajectory has the 
following form [UJ HQ [25] 

7^ = [2T cusp (g) + e(g,j)]lnN + 0(l/L), (1.10) 

where j ' = l/£ = L/(\nN) <C 1 and the scaling function e(g,j) is given by a double series in g 2 
and j. From the point of view of the Bethe Ansatz, the asymptotic behavior (11.101) corresponds 
to the configuration of Bethe roots when all but two intervals [cm-i, cr 2 j] are shrunk into points. 
It proves convenient to describe the asymptotic behavior (jl.lOp using a complimentary set of 
L parameters Sf. (with k = 1, . . . , L) given by the roots of the transfer matrix, the so-called 
'holes'. Then, the relation fll.lOf) corresponds to the configuration [T7] when two holes are large, 
Si ~ — 5l ~ N while the remaining L — 2 'small' holes condense on the interval including the 
origin of length ~ L/ln N. In this way, the factor 2 in front of the cusp anomalous dimension 
on the right-hand side of (11.101) just counts the number of 'large' holes while the scaling function 
e(g,j) describes the dynamics of 'small' holes. 

Examining the excited trajectories close to the ground one, one finds that there exists a 
special subclass of trajectories for which the leading large- N behavior is 7^ ~ (K + 2) IniV 
with K = fixed as L, N — > 00. In this case, the Bethe roots condense on K + 2 intervals while 
the remaining L — K — 1 intervals collapse into points such that [TT1 [26] 

7 JS = [{K + 2)r cusp G7) + e(g,ji, . . .J K+1 )} IniV + 0(1/ L) . (1.11) 

Here ji = a«/£ = EiL /(N In N) are scaling variables (with a« = £i/N being the filling fractions 
of Bethe roots on the cuts) obeying ^f^ 2 ji = L/lniV = j. Similar to (11.101) . the asymptotic 
behavior (II. lip corresponds in the dual, holes description to the configuration in which (K + 2) 
holes are large. For different K the scaling functions e(g,ji, . . . , Ja'+i) are related to each other 
through the Whitham flow which describes a continuous deformation of the distribution density 
of Bethe roots/'small' holes defined for different values of the scaling parameters jk [22j [27] . 

We would like to stress that for large iV the number of trajectories inside the band (II. 6p scales 
as N L ~ 2 and the relation (11.111) only holds in the lower part of the band (II. 6p . The asymptotic 
behavior of the anomalous dimension close to the upper boundary in (11.61) has a different form 

7^ = ^Ws) IniV + 7L (In N/N, l/N) , (1.12) 

where the function 7^ is given by a series in 't Hooft coupling with the perturbative coefficients 
determined in turn by a series in In N/N and l/N. In this paper, we will compute ji for L = 3 
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to three loops in M = 4 SYM theory and will demonstrate that the dependence of 7^ on In N/N 
has an interesting iterative structure. Namely, the perturbative coefficients in front of powers of 
haN/N to higher loops can be expressed in terms of the cusp anomalous dimension [28] 

7^ = £ r cus P G?) In (N + ±LT cusp (g)) +... (1.13) 

where ellipses denote subleading corrections suppressed by powers of 1/N. This property can be 
made explicit by introducing a new scaling function fi{N) related to the anomalous dimension 
via the recursive relation 

7fe x) = h(N+ \in,l) = f L (N) + l(fl(N))' + Ufl(N)Y + ■■.. (1-14) 

Here the prime denotes a derivative with respect to N. We will demonstrate by explicit cal- 
culations that for L = 3 the large- N expansion of the function /i(iV) has a form similar to 
(I1.12p with the only important difference that subleading corrections to fi{N) do not contain 
In N/N l terms (with i > 1) to three loops at least. This implies that In N/N e corrections to the 
anomalous dimension ^^l^ are induced by the leading term in the large N asymptotics of the 
scaling function fi,{N) ~ LT cusp (g) hxN + . . .. The same iterative structure has been previously 
detected for the anomalous dimension of twist two [T9j [29| I3"0l [28] and for the minimal anoma- 
lous dimension of twist three [31] . in which case the corresponding scaling function ceases to be 
In N/N l independent starting from three loops. 

Presently, only the lowest trajectory (11.101) received thorough studies. As it is evident from 
Eq. ( ll.iop . the leading logarithmic asymptotics in the large- N limit allows one to restrict con- 
siderations to the study of the coupling constant dependence of the cusp anomalous dimension. 
The latter was found as a solution to an integral equation for the density of Bethe roots [321 [33] 
or the resolvent associated with the Baxter function [31] . The perturbative solutions to the cusp 
equation were confirmed by field theoretical calculations to four-loop order in 't Hooft coupling 
[2"T] , with the latter providing numerical confirmation for the conjectured all-order integrable 
structures in maximally supersymmetric gauge theory. Many attempts [351 ESI [371 EH EH 00] 
to find a systematic strong coupling solution to the aforementioned equation culminated with 
the development of a scheme which allows one to develop inverse coupling expansion |41j (see 
also [42J). Due to the strong/weak nature of the AdS/CFT correspondence, this provides a 
prediction for the spectrum of corresponding string states. Namely, the anomalous dimension 
(ll.lOP is related to the energy of a classical folded string spinning with large angular momenta 
N and L on the AdS3 x S 1 part of the target space [121 13] • The quantum corrections to the 
string sigma-model were calculated to two-loop order [13] and were found to coincide with gauge 
theory prediction for the cusp anomalous dimension at strong coupling [UJ. In addition, the 
scaling function e(g,j) has been conjectured to be related to the energy density of the ground 
state of bosonic 0(6) sigma-model [25] and the same result was obtained on the gauge theory 
side [4*1"1 |4"5], thus providing nontrivial dynamical test of the conjectured AdS/CFT correspon- 
dence. The string configurations describing the trajectory (11.111) were suggested to correspond 
to a string with K + 2 spikes approaching the boundary of the anti-de Sitter space [131 [TBI SSI I2"B] 
and rotating with large angular momentum N. 

So far, no attempts were made however to address trajectories corresponding to excited states 
interpolating between the low and upper boundaries of the spectrum on the gauge theory side, 
on the one hand, and to find string configurations dual to them, on the other. Thus the goal 
of our present analysis will be to perform the first step in this direction. We will study the fine 
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structure of the spectrum of high-twist operators inside the band in the first few lowest orders 
of perturbative expansion in 't Hooft coupling. Apart from being of interest to the gauge/string 
duality, the considerations done here elucidate an interesting iterative structure of higher twist 
anomalous dimensions. 

Our subsequent presentation is organized as follows. In the next section we outline our formal- 
ism based on the all-order asymptotic Baxter equation in the SL(2) subsector in the maximally 
super symmetric gauge theory. Then we analyze solutions to the Baxter equation corresponding 
to states close to the lower and upper boundaries of the band using two complementary methods, 
asymptotic and quasiclassical considerations, respectively. We find that both techniques allow 
to describe the fine structure of the spectrum for trajectories going deep into the interior of 
the band. They are discussed in Sections [3] and HI respectively. Finally we formulate iterative 
relations for the anomalous dimensions of higher twist operators and then, we conclude. 

2 Asymptotic Baxter equation 

We focus in this paper on the rank one sector of the M — 4 super Yang-Mills theory, which is 
autonomous to all orders in 't Hooft coupling constant. It is spanned by single trace operators 
built from the complex scalar fields X and covariant derivatives 

Otv,l(0)= c fcll ..., fc£ tr{p* 1 A'(0)...Pj i X(0)} . (2.1) 

k!+...+k L =N 

Here the + subscript stands for the light-cone projection of the corresponding Lorentz index, i.e., 
T> + = Dyji^ with n 2 = 0. These operators are transformed into each other under the collinear 
SL(2) conformal transformation of the four-dimensional conformal 50(4,2) group. 
The operators (12. ip obey the renormalization group equation 

+ 7iV ' L ^) ° N ' L ^ = ° ' (2 - 2) 
with anomalous dimensions admitting an infinite loop expansion in coupling constant, 

oo 

7(<7) = J> 2n 7n-i ■ (2.3) 

71=1 

Each term of the series possesses a nontrivial dependence on the quantum numbers L and N 
and, as a consequence of the exact solvability, receives an additional dependence on a set of L — 2 
integers I introduced in Eq. (11.91) . Below we review integrable structures arising first at one loop 
and then extend them to all orders in 't Hooft coupling within the framework of the Baxter 
equation. 

2.1 Leading order Baxter equation 

To start with, let us consider the one-loop anomalous dimensions. Integrability allows one to 
map the dilatation operator D = d/dlnfi into the Hamiltonian of the noncompact spin-^ magnet 
such that the leading order anomalous dimensions 70 coincide with its eigenvalues. To find the 
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latter we resort to the method of the Baxter operator which is known to be naturally adopted 
for analyses of spin chains with large values of quantum numbers. 

The central role in the construction is played by the Baxter Q-operator which depends on 
the complex spectral parameter it, acts on the spin chain sites and obeys a second order finite- 
difference equation. The latter reads for the noncompact spin-| representations living on the 
chain sites (17] 

(u + ) L Q(u + i) + (u-) L Q(u - i) = t(u)Q(u). (2.4) 

Here we introduced a shorthand notation for the imaginary shifts proportional to the value of the 
spins in the spectral parameters = tt± |. The right-hand side of the Baxter equation involves 
a new operator t(u) which is identified as the auxiliary transfer matrix. It is a polynomial of 
degree L in the spectral parameter u with coefficients defined by the local charges q^ acting on 
the spin chain sites and simultaneously commuting with the Hamiltonian, 

t(u) = 2u L + n_ 2 u L - 2 + • ■ ■ + q L . (2.5) 

The definition of the charges q^ is ambiguous since one could have chosen to expand t(u) around 
some reference u = Uq. A distinguished feature of (12.51) is that it does not involve (9(m l_1 ) term, 
or equivalently qi = 0. As we will see later, qi receives perturbative corrections starting from 
two loops. 

The Baxter Q-operator is diagonalized by all eigenstates of the magnet. While the construc- 
tion of spin-chain eigenstates requires the operator itself, finding of the energy spectrum needs 
only the knowledge of its eigenvalues which we denote bjH Qo(u). Since the Baxter operator 
commutes with all local charges, the equation for its eigenvalues takes the form identical to Eq. 
(12.41) . This equation however does not fix the form of the function Qo(u) and has to supple- 
mented by the condition on its analytical properties. In the present setup this boils down to the 
polynomiality of Qq{u) in the spectral parameter u. The large-w asymptotics stemming from the 
Baxter equation fixes its degree d = N to be related to the eigenvalues g 2 ,o of the integral of 
motion q 2 . The latter is related to the quadratic Casimir operator of the s/(2) algebra acting on 
the entire spin chain with eigenvalues related to the total conformal spin of the chain, N + |L, 

g 2 ,o = —(N + \L){N + \L-l)-\L. (2.6) 

Being a polynomial of degree N, Qo(u) can be parameterized by its roots Uk,o as 

N 

Qo(u) = Y[(u-u kfi ) ■ (2.7) 

k=l 

Since Eq. (12.41) is a homogeneous equation, the overall normalization in Eq. (12.71) is inessen- 
tial. The solution to it simultaneously gives the zeroes Uk t o of Qo(u) and quantized values for 
the conserved charges g^o- These allow one to determine the energy, i.e., one-loop anomalous 
dimensions, and the quasimomentum as 

7o .i( ln o4±i)V, * = - iln *M (2 . 8) 

2 V Qo(-5)/ Qo(~k) 



2 We dress the eigenvalue with the subscript to designate the fact that we are considering only leading order 
of the perturbation theory. Above one-loop, the operator formalism which generalizes the one for the short-range 
spin chains based on the existence of ]R-matrices is not yet available. 
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Figure 1: The eigenvalues of the energy 70 and the conserved charge q^o = q3,o/(N + |L) 3 for 
three-site spin chain. For each value of the conformal spin N there are m = ~(iV — 1) + |mod(iV — 
1,3) eigenvalues of the anomalous dimensions with zero quasimomentum 6q = 0. 



respectively. The cyclicity of the single trace operators (12.1 j) imposes a selection rule 9 Q = on 
the quasimomentum. 

Substituting the eigenvalues (12.71) into the Baxter equation and taking the residues of both 
sides at u — Uk,o, one immediately finds that the zeroes of Qo(u) obey a set of N transcendental 
equations 

U kfi \ TT U k,0 — Ujn 



n , (2.9) 

which are identified with Bethe equations. 

The Baxter equation (|2.4j) can be solved numerically for fixed values of the parameters L 
and N. We display in Fig. [1] the result of such numerical solution for three-site spin chains for 
eigenstates with zero quasimomentum and the values of the conformal spin varying in the range 
< N < 50. The following comments are in order regarding these results. For each eigenvalue 
of the energy there are two values of the conserved charge q^^, so that the former are double 
degenerate 

7o(<?3,o) = 7o(-?s,o) , (2-10) 
except for the ground state for which g 3j0 vanishes. A naked eye inspection unintentionally traces 
trajectories in the spectrum in Fig. [lj As it is known from previous studies [221 El IH El E], 
this is indeed the case as eigenvalues fall into families of trajectories enumerated by an integer 
I. For chains with L sites the trajectories are parameterized by a set of L — 2 numbers I = 
(£ 1 ,£ 2 ,...,£l-2). 



2.2 All-order Baxter equation 

Going beyond leading order of perturbation theory one immediately finds that more than two 
nearest-neighbor fields simultaneously interact with each other in the local Wilson operators (12.11) 
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such that the dilatation operator becomes long-ranged. Though, the resulting structures could be 
potentially embedded into higher Hamiltonians stemming from the fundamental transfer matrix 
of the same short-range SL(2) spin chain, this simple possibility does not realize. Therefore, 
one looses operator formulation of the spectral problem based on the existence of the M-matrix. 
Nevertheless, one can consistently deform the leading order Baxter equation to incorporate higher 
loop corrections in 't Hooft coupling. The string-theory considerations imply that the spectral 
parameter gets renormalized beyond leading order of the 't Hooft expansion, i.e., it acquires the 
coupling constant dependence [32], 



x\u\ 



\(u+y/tf=f). (2.11) 



The construction of the deformed Baxter equation relies on the function Q(u) which is still 
polynomial in u but its roots develop a nontrivial ^-dependence, 

N 

Q{u) = Y[(u-u k (g)). (2.12) 

fe=i 

The resulting equation preserves the form of the second order finite-difference equation [4"9l I3"l] 

(x + ) L e A ^ x+) Q(u + + (x~) L e A -^Q(u - i) = t(u)Q(u) , (2.13) 

however with rather involved dressing factors depending on x ± = arfiz*], 

A ± («) = a ± {u) - e(u) . (2.14) 

Here 

- 1 dt lnQ(±i-gt) ( _ ^~^g 



°±(u)= — z 7f= f i- v . ; > (2-15) 

J -i n y/l-t 2 y u + gt 

and Q(u), known as the magnon dressing phase [50] [51], [33] , is responsible for the smooth inter- 
polation between the weak- and strong-coupling expansions [52J, 

KJ "J-iVT^ <?(+§-</*) 7-i s-t 

x f — 1 ln(l + ^ r J ( 1 ^ ) . (2.16) 

Jc [itioo] 2m sinh (7r«) V Axx[K + gs\J\ Axx[n-gs\J 

In complete analogy with the one-loop case, evaluating both sides of the Baxter equation (12.131) 
at u — Uk{g) and taking into account that the renormalized transfer matrix t{u) is regular at 
these points one obtains the long-range Bethe Ansatz equations of Ref. [9]. 

For the fixed length L of the operator (12.11) . the Baxter equation (12.131) fails to encode loop 
corrections above the order g 2 ^ L+l \ When translated to the range of interaction in the dilatation 
operator this happens when the interaction starts to wrap around the chain of fields in the Wilson 
operator (12.11) . This is the so-called wrapping problem which endows the Baxter equation merely 
with asymptotic character. In the lack of the operator framework, the functional form of the 
auxiliary transfer matrix does not appear to be constrained presently by symmetry considerations. 
However, the necessity to match the solutions to the Baxter equation and perturbative results all 
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the way to the wrapping order, i.e., including g 2( - L+1 \ imposes severe constraints on the spectral 
parameter used as the argument of t(u). Namely, t(u) admits the form 



t(u) = Ke 



(x 



(2.17) 



fc>0 



with the sum extending beyond k = L. This implies that the transfer matrix develops non- 
polynomial effects in the spectral parameter. However, these are required to precisely cancel 
the ones emerging from the dressing factors in the left-hand side of Eq. (12.131) . The expansion 
coefficients in Eq. (12.171) are all real, 3mOj;(g) = with Qo(g) = 2. Notice also the appearance 
of the coefficient Hi(g) which was absent at the leading order. 

As we just pointed out, the transfer matrix ( 12.171) is necessarily nonpolynomial in the renor- 
malized x, and after expansion in 't Hooft coupling constant, it preserves this virtue also in the 
bare spectral parameter u. However, if one limits considerations to (and including) order g 2L of 
perturbation theory for the anomalous dimensions fl 2 . 3 j) . one finds that nonpolynomial terms in 
bare u vanish and the transfer matrix can be cast in the polynomial form, 



t(u) 



(2.18) 



fc=0 



with qo(g) = 2. The perturbative expansion of charges in both transfer matrices, Eqs. (12.171) and 

QM = ]>> 2 Xn, Ofc(ff) = ]>> 2n n fc ,n, (2-19) 

n>0 n>0 

allows one to establish order-by-order relations between them for k — 0, . . . , L 

k 

m,0 ) 



Qk,0 
Qk,l 
Qk,2 

with 



m=0 
k 



Ck,m( L ) -\{L-m + 2)^ m _ 2 n m , } , (2.20) 

m=0 
k 

^ C k ,m{L) {0. m ,2 ~ \{L -771 + 2)6 m _ 2 0.m,l + ^{L ~ 777, + 1)(L - 771 + A)6 m ^O. m fl} , 



m=0 



c k,m(L) 



~)k—m 



l-T) cos (f (fc - m; 



(2.21) 



While for k > L, one finds 



Qk>L,0 — 


o, 








= 




Hfc>L+l,0 


= o, 




&L+1,2 = 


\0-L-2,l , 


O-L+2,2 = 


JqO-L-2,0 , 


Q-k>L+2,2 — 



(2.22) 



The charges qi{g) and 92 (p) can be found in a closed form in terms of the Baxter polynomial 
Q(u) to all orders of perturbation theory. Using Eqs. (12.151) and (12.161) . one finds that the first 
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few terms in the 1 / u-expansion of the dressing factors read 

A±(u) = ^± + + 0(u- 3 ) , (2.23) 

where the residues are (a = 1,2) 



u u 2 



Ai a) = v^L^(gt,g)(\nQ(±^-gt))'-^(gt,g) (in ^f^ )'} , (2.24) 



and the explicit form of the functions entering the integrand is 



w^(t,g) = -g\ &)( t ,g)=32it?^, (2.25) 

T 

wW(t,g) = y 2 t, #W(t,g) = 8iZ 2A (g), 

with 

^>-(!)T*^- 

Then one immediately finds by matching the Taylor series of the left- and right-hand sides of 
Eq. (12.1 3p as it — > oo that 

qi (g)=A^(g) + A ( -l\g), (2.27) 

fcfo) = -C 2 (g) + \ {^l\g) + A W ( 9 )) 2 + A?^) + (g) - \{l + 2g 2 )L . (2.28) 

The first term in q 2 (g) is the quadratic Casimir of the sl(2) algebra 

C 2 {g) = (JV + ±L + §7(<?))(Ar + iL + § 7 (#) - 1) , (2.29) 

renormalized by the anomalous dimensions of Wilson operators with the conformal spin iV + ~L, 

1 (g) = -2t^\g)-A W (g)) . (2.30) 
The eigenstates have to be supplemented by the condition of the vanishing quasimomentum 

« = I f 1 * to «±iz4 = o. (2.3i) 



Finally, before we close this section let us point out yet another (very suggestive) form of the 
auxiliary transfer matrix. It allows one to eliminate the superfluous charge q% by renormalization 
of the transfer matrix t(u) via 

t(u) = eh A +^ +A ^ T {u) , (2.32) 

such that 

\ + £ ) , (2.33) 



T (u) = x L 



k>2 



and C\ 2 (g) becomes identical to the eigenvalues of the renormalized quadratic sl(2) Casimir op- 
erator, 

q 2 (g) = -C 2 (g)-\L. (2.34) 
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2.3 Analytic structure of the spectrum 



As we already emphasized in the Introduction, the anomalous dimensions of higher twist opera- 
tors occupy a band of width (L — 2)r cusp (g) IniV for large conformal spins N ^> L. Numerical 
analyses demonstrate remarkable regularity of spectra, see, e.g., Fig. [H of either the energy or 
the conserved charges. Understanding of these properties naturally emerges from semiclassical 
considerations, with the Plank constant played by the parameter 

7 I - 1 = N+\L. (2.35) 

In the small- 77 limit, the Baxter equation reduces to the Shrodinger-like equation for a particle 
in a singular solvable potential. The size of quantum fluctuations is suppressed for vanishing r\ 
and the quantum mechanical motion of particles is restricted to their classical trajectories with 
finite periods. The classical motion governed by the conserved charges represents a nontrivial 
solitonic wave propagating on the chain with L sites. Their quantization stems from imposing the 
Bohr-Sommerfeld conditions on periodic classical trajectories. Thorough studies of leading order 
Baxter equations have demonstrated that the eigenvalues of the Hamiltonian (or the anomalous 
dimension mixing matrix, in the basis of local Wilson operators) fall onto trajectories parameter- 
ized by a set of integers I. This property can be understood within the quasiclassical framework 
as a consequence of the adiabatic deformation of the solitonic (finite-gap) solution with continu- 
ously varying iV which, according to the adiabatic theorem, will belong to the same trajectory. 
Since the Hamiltonian is hermitian, the trajectories do not cross each other as one varies the 
spin N. The hermiticity of the mixing matrix is preserved beyond one-loop order of perturbation 
theory thus one anticipates that these trajectories, while acquiring nontrivial dependence on the 
coupling constant, stay non-intersecting for changing N and g. This implies that the numbers £ 
remain good quantum numbers to any loop order and eventually at strong coupling. 

It is known [22J that the Bethe roots of the short-range SL(2) magnet with L sites, which 
parameterize the zeroes of the Baxter polynomial, admit real values only and populate {L — 1) 
finite interval on the real axis in the quasiclassical limit 77 — 0, 

u E [ax, a 2 ] U [0-3, 04] U • • • U [a 2 L-3, 02L-2} • (2.36) 

The leading asymptotic behavior of trajectories in the small //-limit arises from configurations 
when Bethe roots condense on just two intervals on the real axis. These correspond to the regions 
of allowed classical motion in the separated variables. From the point of view of the spectral 
curve encoding the analytic structure of particle quasimomentum, this corresponds to just two 
open finite cuts and the rest being collapsed into points, i.e., they do not have any Bethe roots. 
The shrunk cuts were dubbed double points in Ref. [T7|. The fine structure of the band emerges 
when these double points open up into cuts. The single logarithmic asymptotics in the large-spin 
limit arises when these two cuts collide at the origin and may or may not pinch the double points 
between them. The low boundary of the spectrum corresponds in this picture to the case when 
all double points are trapped at the origin by the inner end-points of the colliding cuts. While 
the upper boundary of the spectrum is a consequence of the migration of all these double points 
outside of the merged single cut involving the origin. The intermediate states in the spectrum 
are described by continuous deformation of the multicut configuration, when the external cuts 
gets shrunk with internal cuts opening up at the same time. 

As our discussion suggests, the upper and lower boundaries of the spectrum corresponds to 
different analytic structure of the distribution of Bethe roots in the complex plane. These can 
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be translated into the asymptotic behavior of the conserved charges qk{g) or, equivalently, the 
roots Sk(g) of the transfer matrix, 

L 

t(u) = 2Y[(u-6 k (g)) , (2.37) 

fe=i 

which similarly to qk{g) admit perturbative expansion in 't Hooft coupling constant 

5 k (g) = J29 2n h,m (2.38) 

ra=0 

and are in one-to-one correspondence to each other, 

L L L 

I>fo) = -§91(0) , J2 5 l(9) = U(9)- 12(9), ... \{h{g) = \{-i) L qL{g). (2.39) 

k=l k=l k=l 

The upper boundary of the spectrum corresponds to large values of all conserved charges, 

9M=V~%{9), (2-40) 

as r] — > and all rescaled charges qk(g) being of order one, qu{g) ~ 0{rp). The same is applicable 
to the roots Sk{g). Thus all quantum numbers determining the spin-chain state are large and 
quasiclassical approximation is well defined in this regime. 

On the other hand, the lower boundary of the spectrum (jl.6p corresponds to the degener- 
ate situation when all higher conserved charges qk>2{g) take anomalously small values, as was 
demonstrated in Ref. [T7| at leading order of perturbation theory, 

Q2k ~ 2^|> + Otf) , (2.41) 

where (5 = |L/ (N + |LJ vanishes in the single-logarithmic limit N ^> L. This indicates the 
appearance of two types of roots of the transfer matrix, 'large' and 'small' possessing the following 
asymptotics in the leading order in 't Hooft coupling [T7j 

\5i,l\ >j= 2 N + 0(g 2 ) , < ^ + 0(g 2 ) , (2.42) 

which exhibits the emergence of a new "hidden" parameter 

£ = -jrhiiV. (2.43) 

For the lower part of the spectrum, the semiclassical analysis is not applicable. One has to reply 
therefore on a complementary method which we will employ below. It will be applicable for 
the spectral parameter u ~ 0(r] ) thus yielding a valid expression for the anomalous dimensions 
which can be used for both the lower and upper parts of the spectrum. 

Our subsequent consideration will be focused on developing different techniques to analyze 
the excited states in the vicinity of the lower and upper boundary of the band. The consideration 
will be limited to the lowest three orders of perturbation theory, though it can be generalized 
in a straightforward fashion to even higher loops. At the lower boundary, the analysis will be 
performed in great generality for any N and L, while for the upper part we will consider the 
three-site spin chain which, on the one hand, is simple and, on the other, demonstrates all of the 
salient features of the method and emerging integrable structures in higher orders of pertubration 
theory. 
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3 Lower part of the spectrum 



Let us begin with the description of the lower part of the spectrum in the limit of large spin 
N. The asymptotic solution to the Baxter equation in this region will be valid for the spectral 
parameter behaving as u ~ O(rj ). The method relies on an observation that two terms in the 
left-hand side of the Baxter equation (12.131) have different scaling behavior as rj — > 0. This method 
has been developed in Ref. [53] for short-range spin chains and applied for the calculation of the 
one-loop minimal trajectory in Ref. [T7]. Our main goal is to extend this approach to higher 
loops. Let us recapitulate the essential features of this formalism at one loop and then extend it 
to arbitrary orders of perturbation theory. 



3.1 Asymptotic solution of one-loop Baxter equation 

We notice that the evaluation of the leading order anomalous dimension (I2.8P requires the knowl- 
edge of the Baxter polynomial and its first derivative at the points u — ±|i. In their vicinity, 
due to the power suppression by the dressing factors, one can neglect either the first or second 
term of the one-loop Baxter equation. This however persists in a more general situation when 
u ~ 0(i] ). In the asymptotic regime under consideration, the transfer matrix is large \t(u)\ ^> 1 
since the conformal Casimir scales as a second power of r/ _1 . Then, introducing the ratio of the 
Baxter polynomials <f(u) = Qq{u + i)/Q (u), one realizes that in order to match the scaling be- 
havior of the left-hand side it should be either large or small 0. Due to this fact the second-order 
finite difference Baxter equation (12.41) at one-loop level splits into two first-order finite-different 
equations [TF] 

( M ± ) L g ± ( M ±0 = t ( M )g ± ( M ), (3.i) 

whose solutions have the form (up to an overall normalization factor), 

Qt(u) = 2^ U ^ T{TiU±i * k > o) , (3.2) 

r(=F*i* + !) 

where we used the parameterization of the transfer matrix at leading order in terms of its roots 

L 

t (u) = 2l[(u - 5 kfi ) . (3.3) 

k=l 

The asymptotic solution of the one-loop Baxter equation is given by a linear combination of 

QiT\u) = A+Q+(u) + Ao Qo («) , (3.4) 

where are fixed up to an overall normalization from the condition of the cyclic symmetry of 
the eigenstates 



exp(^ ) = = 1 =► i4± = Qf( T i<) ■ (3-5) 

vo(- 



3 For the spectral parameter scaling as u ~ 0{r\ 1 ) both terms in the left-hand side would contribute on equal 
footing. 
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3.2 Asymptotic solution of all-loop Baxter equation 

The asymptotic solution to the one-loop Baxter equation can be generalized to higher loops. To 
this end we have to assume that higher loop corrections to the conserved charges qk{g), roots 
of transfer matrix Sk(g) and Bethe roots are compatible with the large- N scaling of the 

corresponding one- loop quantities, i.e., higher loop corrections should not violate their scaling 
behavior at one loop so that \t(u)\ ^> 1 order-by-order in 't Hooft coupling. Analogously to the 
leading order equation, expanding both sides of the Baxter equation ( 12. 13ft in powers of 't Hooft 
coupling, one observes that in the vicinity of the points u = ±~z one can neglect either set of 
terms arising from the expansion of the two terms in the left-hand side of the all-loop Baxter 
equation (12. 13[) . A close inspection shows that, contrary to the one-loop case where it is valid 
for any length L of the chain, this assumption is fulfilled starting from certain values of L which 
is larger at higher orders of perturbation theory. 

Having these limitations in mind, the asymptotic solution to the all-loop Baxter equation 
( 12. 13[) can be split into two equations 



(x ± ) L e A±(u±) g ± (w ± i) = t{u)Q ± {u) . (3.6) 

Now we proceed with the solution of these equations in the first three orders of the loop expansion. 
Taking into account that the only dependence of the Baxter polynomial and the transfer matrix 
on the coupling constant comes through their roots, the perturbative expansion of the latter 
reads 

Q(u) = 9 2n Qn(u) , t(u) = 9 2n tn(u) , (3.7) 

n=0 n=0 

where the leading order transfer matrix t (u) is given by Eq. (13. 3p while the higher terms are 



ti(«) = - «o(«0. (3-8) 




j<n U ~ 6 i'° U ~ " n '° / 

Substituting these relations into the asymptotic Baxter equations (13.61) and expanding the dress- 
ing factors to the required order, we deduce the first-order finite difference equations for Qi(u) 

(u^Qfiu ± i) - ( M ± ) i - 2 (\L =f m ± 7o ± ) Qo(u ±i) = Qf(u)t (u) + Q^{u)h{u) , (3.9) 

and Q 2 (u) 

(u^Qfiu ± i) - (^) L - 2 {\L T iu^t) Qf(u ± i) (3.10) 
+( M ± ) i ~ 4 (l 2 L{L — 3) =f |(L - 1)^ T i^fof ± nti^) 3 ) Qo( u ± 

= Qt(u)t (u) + Qf{u)h{u) + Qt(u)t 2 (u) . 

In these equations 7^ stand for the coefficients in the perturbative expansion of the anomalous 
dimensions (I2.30|) 

7(s) = $> 2n+ V, 7, = 7„ + -7", (3.H) 

n=0 
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with 7^ associated with A±\g) in (12.301) . and the constants af are defined as 

3i ( Qo(±f)\ 2 i Qo(±|) 



1 i6 vg (±|)y i6Qo(±|) 1 ; 

As can be seen from the structure of the asymptotic Baxter equations (13.91) and (13.101) . beyond 
one loop level it is natural to look for the solution for Qf (u) in the following factorized form 

Qf(u)=Qf(u)Gf(u), n>0. (3.13) 

The solutions are easily constructed in term of Euler psi-functions and read, 

L 

Gf(u) = cf + \LiP'{Tiu + §) - Jo^iTiu + f) ± i J2 <W(=R« ± i6 k)0 ) , (3.14) 

k=l 

L L 

Gf(u) = c± + § (Gf(u)) 2 - \ ± %o) ± i 6 iMTiu ± % ) 

3=1 3=1 

- ±btf"(Tiu + \) + + |) + (±iaf - |( 7o ± ) 2 ) ^(T«« + \) ~ 7^(T^ + §) . 

Note that the two and three- loop asymptotic functions G l2 {u) contain arbitrary constants c 12 
reflecting the fact that the solutions of the asymptotic Baxter equation (13. 6p are defined modulo 
multiplication by an arbitrary constant Q ± (u) — > c ± (g)Q ± (u) with c^{g) = 1 + cfg 2 + cfg 4 + . . .. 
As we will see momentarily, although the constants cf 2 appear in the solutions, they do not 
contribute to the anomalous dimensions. 

The asymptotic solution to the Baxter equation is given by a linear combinations of Q ± {u) 

Q^(u)=A + (g)Q + (u)+A-(g)Q-(u), (3.15) 

with the coefficients A ± (g) admitting perturbative expansion in the coupling constant 

A ± {u) = Af(u) + g 2 Af(u) + g 4 Af(u) + 0(g 6 ) , (3.16) 

and with analogous expansions for Q^(u) 

Q (u) = A+Q+(u) + AqQq(u) , (3.17) 
Qi{u) = A+Q+(u) + A^Q (u) + A+Qt(u) + A Q^(u) , .... 

Recall that to one loop, the normalization coefficients Af are fixed by the cyclicity condition 
exp(z0) = 1, Eq. (I3.5p . In a similar manner, constraints on the normalization parameters Af 2 
come from the fact that the quasimomentum 9 is protected from perturbative corrections and, 
therefore, it does not depend on the 't Hooft coupling constant g. 

3.3 Asymptotic anomalous dimensions 

Using the expression for the anomalous dimensions in terms of the all-order Baxter polynomial 
(I2.30p . which explicitly reads 



7 (g) = ig 2 £ f VT^T* (in %|+f _ ^ ) =7 + (s)-7-(s) = 2Ke 7 + (s), (3.18) 
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one immediately finds the first three terms in 't Hooft expansion 

7o ± = |(lnQ (±f)) / , (3.19) 
7i ± = |(Gf(±f)) , + ^(lng±(±f)) /,/ , 

7^ = | (G±(±i))' - |C±(±i) (G±(±§))' + ± (Gf(±l ))'" + gij (lng±(±| )) (5) . 

Substituting the asymptotic solutions (13.21) and (I3.14p into (13.191) we immediately arrive at the 
following expressions for the anomalous dimension in terms of the roots of the transfer matrix S k 

L 

70 = In 2 - 14(1) + Ke^VK! + *M , (3-20) 

k=i 

L L 

71 = |L<(l)-|7o^(l)-|^e^V"(i + z4,o)- Sm^4,/(i + *,o) , (3.21) 

fc=i fc=i 

72 = |^(l)7o - |^(l)7i - iSfe (( 7o + ) 2 - 2*a+) V"(l) - |^ {4) (1) + ik^E^I + W *.o) 

fc=i 

L L L 

+ §3m ]T + iSk,o) ~ |3Re ]T ^"(f + z5 fe)0 ) -3m^ 4, 2 ^'(| + iS k>0 ) . (3.22) 

fe=i fc=i fc=i 

We see that though the solutions to the Baxter equation are plagued by the arbitrary constants 
cf 2 they disappear from the expressions for the anomalous dimensions (I3.20p - (I3.22p . For the 
lowest trajectory in the spectrum with q 3 = 0, we immediately reproduce findings of Ref. [54J. 

3.4 Quantization conditions 

The anomalous dimensions (I3.2(jp - (I3.22p depend on the roots of the transfer matrix. In this 
section we will find their quantized values which yield in turn the quantized charges q% parame- 
terized by L — 2 quantum numbers £ = {£%, £2, ■ ■ ■ , £1-2}- Each set of the charges gives rise to a 
certain state in the spectrum of anomalous dimensions at fixed spin N and twist L. 

In the method of the Baxter Q-operator the quantization conditions follow from the require- 
ment for Q(u) to be a polynomial in u. In the large- N limit, or equivalently for r) <C 1, the 
Bethe roots scale as ~ N ^> 1 and therefore we are not allowed to impose the polynomial- 
ity condition on Q^ s \u) since this solution is only valid for the spectral parameter u ~ 0(r] ). 
Instead, we have to demand that the asymptotic solution Q^^u) has to be a regular function 
on the real axis for u ~ 0(ff). As it follows from (13.131) and (I3.14p . possible poles in Q^(u) 
originate only from the ^-functions on the right-hand side of (I3.14p and their position is deter- 
mined by roots of the one- loop transfer matrix 5kfi- In the lower part of the spectrum, all but 
two roots scale as 5k,o ~ 0(i] ) (with k — 2, . . . , L — 1) while the two remaining ('large') roots 
scale as #1,0 ~ — 5l,o ~ ^(v 1 )- We require that Q( as )(u) should have zero residues at 'small' 
roots u = 5k,o ~ 0(rf), k — 2, . . . , L — 1 and obtain the following quantization conditions 

res {A + (g)Q + (u)+A-(g)Q-(u)} = 0, (3.23) 



u=8, 



k.o 



where the constants A ± (g) are determined from the cyclicity condition on the single-trace Wilson 
operators. 
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At one-loop level, we take into account the relations (13.21) and (13.51) and obtain from ( 13.231) 
the quantization conditions on the roots of the one-loop transfer matrix [TT] 



n 

Mi 



3,0) 



r(tf. 



1(5 



fc,dJ 



r(| + i6 jt o) 1 L -q T(i - i<5 fc , 



r(| -i6j,o)i 



1 . 



(3.24) 



Note that the product over j involves all roots, large and small, while the conditions are imposed 
on the small roots only. For the minimal energy we have only two large roots 8\fi = —5l,o, while 
the rest of the roots are small and paired 5k,o = — ^L-fc+i.o with k — 2, . . . , L — 1. Separating the 
contribution of the large roots and taking the logarithm of (I3.24p we get 



L-l 



5 nfi ln(-g 2 ,o) + £argr(| - iS nfi ) + ^ argT(l + i<5 ni0 - i6 jfi ) 

3=2 



2 " ""n i 



(3.25) 



which can be considered as an equation for the small leading order roots (5 n> Q. The integer 
numbers k<i > k 3 > . . . count the branches of the logarithms and play the role of quantum 
numbers in the system of L — 2 quantized charges. They satisfy the condition k n = —kL- n +i- 
The explicit dependence of integers k n on their index n for the minimal energy will be obtained 
in the next subsection. Expanding Eq. (13.251) in the limit |<5 nj o| < 1, we finally arrive at the 
following quantization condition [TT] 



^nk r 



In(-g 2 , ) + (L — 2)^(1) - 



(3.26) 



The exact two- and three-loop quantization conditions are more involved. To save space, we 
present here only the two-loop result 

3fJe[^ 1 [ln2 - LV(§ - iS j>0 )] - ^ [O^M ~ 5 3,M i5 J,o ~ *4,o) - 4,i^(| - *4,o)] } (3.27) 

= 3m j 7o + (V(l) - - i&jfi)) + W(\ ~ iS J,o) - I E ^(i 
^ k=i 



iSk,o) 



Similarly to Eq. ( 13.261) . these quantization conditions simplify for small higher order roots. Keep- 
ing only the leading terms in their Taylor expansion, they read at two- and three-loop orders, 
respectively, 



21n(V2<5 1 , ) - + {L~ 2)^(1)1 8 n ,i = J (2 7 o^(|) - *n,o + 0{8 2 



21n(V2 ( 5 1)0 )-L^(|) + (L-2)V(l) 



->n,2 



| (2 7o ^(i) - L^'(i)) <S, 



(3.28) 



L0(4)(I) _ 8^(|) 7o 2 - 4^"(|) 7 o + 32 7l ^(i) - 2( X + + x-)^(l) 



M,0 



+U4 - cr)(2 7 o^(l) - ^'(1) + 47r5 ltl - c+ - q) + 0(S 2 ), 



up to higher order terms in small roots, denoted cumulatively by 0(S 2 ). These linear equations 
relate small roots of the transfer matrix at successive orders of perturbation theory with coef- 
ficients accompanying them which scale logarithmically with iV stemming from the anomalous 
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dimensions and logarithms of the large roots of the transfer matrix. To simplify notations, we 
introduced the following combination 



x± . 16(70 ^ ±16ia -(|g|) 2 + g||. (3 , 9 ) 

We notice that at large N the coefficients in front of 8 n> i and 5 n> o in both sides of the first 
relation in f!3.28|) are of the same order O(lniV) and, therefore, £ n) i/<5 n) o = O(N ). We also 
observe that the two-loop constants cf do not show up in the two-loop condition, however, they 
enter into the three- loop quantization condition, the second relation in (13.281) . A quick inspection 
shows that while all c^-independent terms are of order on^|, the contribution proportional to 
cf — c{ scales as In N, i.e., 

2 7o ^(1) - \tx 2 L + 4^1,1 - q - cf ~ 2(2^(1) + V2n) In N > 5 nfl . (3.30) 

We recall that one of the starting points of the asymptotic method was the assumption that the 
large- N scaling of higher loop corrections to the 5-roots should be compatible with the scaling of 
the leading order roots 8k,o- I n order for the three- loop quantization conditions to be consistent 
with this assumption, or equivalently the corrections 8 n ^ to be "small", we have to impose the 
following constraint on the two-loop constants 

cf = ci . (3.31) 

Only in this case our initial assumption validating the current approach will be satisfied. Fixing 
the relation between cf in this way and using in the leading order 

7o ~21niV, 7l ~-±7r 2 lniV, X + + X~ - ~H » ( 3 - 32 ) 
we finally get the two- and three-loop corrections to S n 

6n,i - \n 2 K,o , &n.2 ^ -^7r 4 ^n,o • (3.33) 

Equations (I3.26P and f!3.33|) together with the asymptotic formulae for the anomalous dimensions 
f )3.20p - (13.221) completely determine the spectrum at its lower boundary up to three loops. 



3.5 Ground state and daughter trajectories 

The expressions for the anomalous dimensions (13.201) - (I3.22p in terms of the roots of transfer 
matrices are valid for N, L ^> 1 at the bottom of the spectrum (11.61) . Below we will analyze 
in detail the asymptotic limit £ ^> 1 and find preasymptotic corrections in l/£ to the leading 
logarithmic behavior (11.61) . To derive the explicit dependence of the anomalous dimensions on 
the hidden parameter £, we heavily rely on known scaling properties of the roots of the transfer 
matrix as a function of the conformal spin found in the previous section. We remind that the 
large roots are related to the eigenvalues of the conserved charge ^(fiO, 

61(g) = ~h{g) = (~ l M9)) 1/2 + ■ ■ ■ , (3.34) 

4 The apparent dependence on In 2 N rather than being linear in IniV in the second term in the right-hand side 
of Eq. (|3.28p cancels between 7q and x + + X" contributions, see Eq. (|3.32p . 



19 



where 92 (p) is defined in Eq. (12.281) such that 

6i, = -= + ..., = lniV + ... , 5 1>2 = -±n 2 \nN + ... . (3.35) 

v2 

On the other hand, the small roots of the transfer matrix scale as 5k ~ l/£. Selecting the large 
roots contributions and expanding the rest in terms of \S n fi\ <l,we get for the one- loop minimal 
anomalous dimension the well-known result [lTj 

L-1 

70 = 21niV-2Lln2 + 7C(3)^(<5 M ) 2 + ••• , (3.36) 
and analogously, applying f!3.33j) to the two- and three-loop anomalous dimensions we find 

L-1 



_ y i n N _ | C(3 ) + L ( C(3 ) + y i n 2 ) - (f C (5) - §C(3)tt 2 ) £ (5 k , ) 2 + . . . , (3.37) 

3=2 

^vrMniV + f C(5) + |C( 3 )^ " ^(f C(5) + ^C(3)tt 2 + ^7r 4 ln2) (3.38) 

+(Tc(7) - W c(5)tt 2 - ic(3)tt 4 ) f>£ >) 2 + . . . , 



fc=2 



where the small roots 5^,0 are determined by the quantization conditions (13.261) . The ellipses 
stand for subleading terms with higher powers of 8's. The quantized values of the small roots Sk,o 
depend on the integer numbers k n = —k n ^L+i defined in Eq. (I3.25p . with the ones corresponding 
to the minimal anomalous dimension being k n = L + 1 — 2n, yielding 

j=2 n=2 

The set of relations (I3.36P - (I3.39|) defines the three-loop minimal anomalous dimension for 
£ 3> 1 with the first subleading corrections taken into account. The present formalism allows us 
to calculate further terms in inverse powers of £, by expanding the Eqs. (13.201) - (I3.22p to higher 
orders in small roots 5. We will not perform this calculation further. The anomalous dimensions 
of excited, or daughter, trajectories are determined by the same expressions however with another 
set of integer numbers k n . The first daughter trajectory possesses the same "occupation numbers" 
as the minimal anomalous dimension except for k 2 = — &l-i = L — 1. This state is separated 
from the minimal one by the following gaps 

J IT 2 

A 7o = |C(3)-^-, (3.40) 



ln^iV 



Ltt 2 



^ = -(fC(5)-aC(3)7r ;inajv 
^72 = (^C(7) - WC(5)vr 2 - ^(3)^) Ln 



\n 2 N 



at one, two and three loops, respectively. 
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Substituting (13.391) into ( 13.361) . (13. 37ft and (13.381) we find that the subleading corrections 
to the ground state trajectory develop a nontrivial dependence on a new scaling parameter 
j = £- 1 = L/\nN 

^(g) = [2T cnsp (g) + e(g,j)] In N+ . . . . (3.41) 

As was shown in Refs. [T7J [2U [25], this scaling behavior holds both at weak and at strong 
coupling. In the large-£ limit, which we addressed above, e admits a double Taylor expansion in 
g 2 and j, 

e{g,j) = g 2 e (j) + g^j) + g 6 s 2 (j) + . . . , (3.42) 
with perturbative coefficients taking the form of the power series in j, 

e (j) = -2 ln(2)j + ^vr 2 C(3)j 3 + . . . , (3.43) 
eiU) = [C(3) + K In 2] j + [-|vr 2 C(5) + |A(3)] f + • • • , 

S2(j) = - [f C(5) + ^vr 2 C(3) + i-yin2] j + [f|7r 2 C(7) - ifvr 4 C(5) - gfeAP)] f + • • • 
shown to order j 4 . This result agrees with the findings of Ref. [Hj (see also Refs. |45l[55l l56 | [57]). 

4 Upper part of the spectrum 

In the previous section, we described the lower part of the spectrum close to the minimal tra- 
jectory. Let us turn to the study of the loop effects in the upper part of the band (jl.6p . As we 
discussed earlier, in the asymptotic regime in question all roots of the one-loop transfer matrix 
are large and scale like Sk,o = O(N) implying that conserved charges behave as qk,o ~ V~ h i such 
that semiclassical expansion is legitimate [17J. As in the previous section, we shall assume that 
the scaling of higher loop corrections to the (5-roots is compatible with the one of leading order 
roots. 

4.1 Semiclassical Baxter equation 

In the considered asymptotic regime it is instructive to rescale the spectral parameter 

u = v^ x u , (4.1) 
and introduce the Hamilton- Jacobi "action" function S(u) as follows 

Q(u) = exp (-S{u)) , (4.2) 



n 



where 

N 



S(u)=7) y %2hy(u-riUk(g)). (4-3) 

fe=i 

One assumes that the rescaled charges qk(g) = rj k qk(g) and the transfer matrix 

L 

t(u) = f] L t(u) = 2u L + Y, u L ~ n q n (g) (4.4) 



n=l 
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admit a regular expansion in powers of r\ 

US) = t(9) + V$ ] (9) + 0( V 2 ) , r(u) = rM(2) + V r [1] (u) + 0{r, 2 ) , (4.5) 

with each term having in turn a perturbative expansion in the coupling constant. In the semi- 
classical approach one also assumes that the "action" S(u) is a series in 77 

S(u) = S [0] (u) + V S [1] (u) + 0(r] 2 ) , 

which is convergent and each term in it is uniformly bounded. As was demonstrated in Ref. 
[T7] this assumption is justified at leading order of the perturbative expansion anywhere in the 
band provided that £ < 1 or equivalently IniV < L. In the case £ S> 1, or IniV ^> L such an 
assumption breaks down in the low-energy part of the spectrum and is valid therefore only at the 
upper boundary. When going beyond one loop we assume similar behavior for the power series 
in 77 of higher order corrections in 't Hooft coupling to the action Sn\u). 

Substituting (14.21) and (14.51) into the Baxter equation (12. 13[) and expanding it in the double 
power series in g 2 and 77 we get at leading order in 't Hooft coupling 

2cos(4° ]/ (n))=4P> ( 4 -6) 
-2S£ ] '(u) sin (4 0] '(£)) = + Sf\u) cos (s [ ° ] '(u)) + | sin . (4.7) 

Here r^{u)/u L plays the role of an effective potential. Analogously, at two loops we obtain 

r\ 0] (u) 



2M 0J/ (^)sin(^)j=^p, (4* 



-2S?\u) sin (S^(u)) = - S?'(u)Sr(u) sin (4.9) 

+ (sW'W + 2S?'(u) S W(u) - + cos (s?'(u)) . 

The derivation of subleading WKB corrections at higher loop orders is straightforward but the 
resulting equations are rather cumbersome to be displayed here. 

4.2 Bohr-Sommerfeld quantization 

The Bethe roots of the SL(2) magnet take real values only. Numerical solution of the multiloop 
Baxter equation suggests that this property persists in higher orders of perturbation theory and 
eventually to all loops. In the semi classical limit, 77 — > the Bethe roots condense on finite 
intervals on the real axis [581 [22] so that their normalized distribution density defined as 

dup{u)=7 1 N 1 (4.10) 



vanishes outside of the domain 

S = [cti,<t 2 ] U [£3,(74] U ••• U [0^-3,0^-2] , (4.11) 
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with ijfc obeying the condition a^ L r(ak) = ±2 [22]. Here the interval boundaries ordered as 
<7i < (72 < • • • < oil-2 and one of the intervals contains the origin. The function S'(u) introduced 
in the previous section is related to the root density via 

S'(u)= I dv^L. 



s u ~ v 



It is a double- valued function on the complex w-plane with square-root branching points 

Let us denote by the number of roots at the fc-th interval [a^fc-i, a 2 k\ with k = 1, L—l. 
These integers do not depend on the coupling constant simply because the number of roots at a 
given interval does not change when perturbative corrections are turned on, 



L-l 



J2 n k = N > 0<m,n 2 , ...,n L ^<N. (4.12) 



k=l 



By definition from (14.101) and (14.111) we have for each interval from S 



°2fc 



dupiu) = r]n k . (4-13) 



°"2fe-l 



Expanding the density of roots in ^-series p(u) = p^(u) +r]p^(u) + . . . we obtain for small values 
of nfc = O(N ) in each order of the expansion 



&2k r&2k r&2k 



dup [0] (u) = 0i / dup [l \u) = Tjn k , / dup [j>1] (u) = 0. (4.14) 



°2fc-l J < J 2k-l J a 2k-\ 

These can be rewritten as contour integrals of the resolvents over the cuts [cx^-i, &2k\ on the real 
axis in the complex w-plane, 

?2fc f du 

dup{u) = d> S'{u) , (4.15) 

Jaj 27TZ 

where the contours «j encircle the intervals [<72fc-i, &2k] i n the anticlockwise direction. As a result 
we obtain the Bohr-Sommerfeld quantization conditions at all orders of perturbation theory 

du 

a 2m S ' { ^ = Vnj ■ (4 ' 16) 

The set of integers n = . . . , nL_i} obeying (I4.12p plays the role of the quantum numbers de- 
termining the quantized values of the conserved charges and, as a consequence, the anomalous 
dimensions of operators. These numbers enumerate the trajectories which pass through anoma- 
lous dimensions for integer values of the Lorentz spin iV and parameterize them starting from the 
highest one residing at the upper boundary of the spectrum. These trajectories represent a legiti- 
mate analytic continuation of the anomalous dimensions. Another analytic continuation reflected 
in a yet another set of trajectories encoded in the set of quantum numbers £ = {£%, £%,... , £1-2} 
was analyzed in the previous section. In fact the two sets of integers are related to each other by 
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linear relations, e.g., i = [N/2] — n for L = 3. Since the right-hand side of Eq. (14.161) does not 
depend on 't Hooft coupling g, we immediately conclude that 

These relations hold to any order of the //-expansion and allow us to compute the quantized 
charges qf 1 for given k and j. In the next subsection we will examine the two- and three- loop 
quantization conditions (I4.17P order-by-order in powers of r\. 



4.3 Quantization for three-site chain 

Let us resolve the quantization conditions for the high-spin Wilson operators of twist three, i.e., 
L = 3 and N > 1. In this case, the Bethe roots condense on two intervals. We assume that 
n = O(N ) roots reside on [cti,^] and the remaining N — n roots are on the interval [S^o^]. 
Then in the leading order of the semiclassical expansion, one finds that 

^T(3) = 0, (4-18) 
«i 2m 

which implies that the interval [cri,^] shrinks into a double point, according to the terminology 
of Section 12.31 The double point a* is a solution to the equations 

rf(a*) = 0, \rf\a*)\ = 2, (4.19) 
and defines the critical values of the conserved charge qf^ which takes the values 

3S?o = ±4f • ( 4 - 2 °) 



The two possible signs of the charge lead to the double degeneracy of the energy spectrum 
which is not sensitive to the interchange q^o — > —93,0- We will see below that this degeneracy is 
preserved in higher loops as well. For our subsequent consideration we choose therefore just one 
value qfl = 1/^/27. 

Let us examine the Bohr-Sommerfeld quantization condition (I4.17P involving the first sub- 
leading semiclassical correction Sq (u) to the Hamilton- Jacobi function. Since the integration 
contour a\ encircles a vanishing interval [ai,a 2 ], the function should possess singularities on 
it to yield a nonzero right-hand side in Eq. (I4.17p . This can be demonstrated by opening up the 
double point a* = 1/ ^/l2 into a small interval parameterized by a small parameter e < 1 and 
populated with n Bethe roots 

pi,^ 2 ], °i,2 = °*T\e. (4.21) 

Using the condition defining the end-points of the cuts Tq (S^) = 2, with the interval param- 
eterized by the variable z such that u = a* + ez, we find in the leading order of the Taylor 
expansion in e that 

<g\^^=-^e\ (4.22) 
3,0 v 7 ^ 2 v ; 
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Ql,l o.o 




Figure 2: The exact eigenvalues of charge q\ at two-loop order of perturbation theory. 



Substituting these results into Eqs. (14. 6[) - (j4.7j) . we can find the quantized values of the first 
subleading corrections to g 3i0 from f)4.17p . 



n . 



(4.23) 



The contour ot\ encircling the interval \d\, a%] can be deformed away from it, such that the integral 



is given by the residue of S^'id* + ez) at z = oo keeping the product ez fixed. Extending these 
considerations to higher orders in the ^-expansion, we get semiclassical corrections to the charge 
§3 at one loop in 't Hooft coupling, 



93,0 



- v- 



27 

n + 1] 



(4.24) 



+ r] 



V3 | 
108 ' 



a/3 

(24n 2 + 60n + 43) - rf /Q - 3 1 u ~ 
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8n 3 + A8n 2 + 47n + 32) + 



The procedure can be generalized to an arbitrary high order in rj. 

Turning to higher orders of perturbation theory, one observes that starting from two loops 
the transfer matrix (14.41) develops a nontrivial coefficient in front of the [L — l)-st power of the 
spectral parameter, i.e., the charge q\{g)- The latter however is completely determined by the 
Baxter polynomial (I2.27P such that the n-th perturbative term in qi is defined by (n — l)-st 
order of Q{u). Using the asymptotic solutions of the leading and next-to leading order Baxter 
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equations ( 13.21) and (13.141) . it is straightforward to obtain qi(g) up to three loops, 
1 L 

i=i 

?i,2 = -^9i,i (4.26) 

L L 

+ ^ E (^(1 - ***.«>) - <(f + <ff *,o)) - ^ E - ***,<>) + v>'(f + <w) > 

fc=l fe=l 

where 

4,o = ^ M = J (? + + ^(^)) , (4-27) 

are the leading order roots of the transfer matrix. They are exactly calculable from the leading 
order charges qk$ for any fixed L. Particularly, in the considered WKB regime one can explicitly 
show from Eq. (14. 24ft that all roots are "large" and scale as 5k,o = O(N). For the maximal 
energy and g3 5 o > 0, two of the roots are positive and one is negative, o 10 > 0, d 20 > 0, 5 3 ' < 0. 
Expanding the expression for two- and three-loop corrections q\ >n ( 14. 25ft in the limit rj <C 1, or 

— Tq"! 

|5[ /?7| ^> 1 and using the asymptotic formula for the Euler digamma function 



I °° B 

ib(x) ~ lnx x^ 1 — — —x~ 2h , x ^> 1 



2 ^ 2A; 



fc=i 



we find (for g 3 ,o > 0) 



TV 2 



9i,i = , 9i,2 = • (4-28) 



Notice that q\ changes its sign as q 3 — > —q 3 such that all explicit expressions for q 1 we found 
above are accompanied by the signature factor sign(g 3 ) to accommodate both values of q±. It is 
interesting to note that q\ does not receive any higher order power corrections in rj. Indeed, this 
result has been confirmed by numerical calculations of the charge q\ t x (see Fig. [2]) 

The resolution of the quantization conditions (I4.17P at higher loops goes along the same lines 
as the one-loop calculation and yields the quantized values of corrections to the conserved charge 
q 3 at every order of the semiclassical expansion 



93,1 



+ r/ 3 || (2 7 S 0] (24n 2 + 60n + 43) - 144 7o 1] (n + 1) + 72 7 J 2] + 216n + 9^) + . . . , 

93,2 = ^ ( 4 ^! 0] - t ) " v2 m ( n2 ~ 2 ^ ( 27i ° ] + t) + 487)01 {n + 1} 

-24 T W-12( T P) 2 -^)+.... 

Here 7j ^ is the p-th coefficient in the semiclassical expansion of the k-th order anomalous dimen- 
sion in 't Hooft coupling g 2 , 

l k = lt ] +Vli 1] +V 2 lt ] + ---- (4-30) 
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10 20 30 40 50 !0 20 30 40 50 10 20 30 40 50 



N N N 

Figure 3: The exact eigenvalues of the conserved charge % at one-, two- and three-loop order (left 
to right) of perturbation theory for twist-three operators shown together with a set of selected 
semiclassical trajectories computed from Eqs. (14.24)) and (I4.29p . respectively, for n = 0, 1,4, 7. 



Several comments are in order. 

We observe that the semiclassical expansion of the two- and three-loop corrections to charges, 
i.e. qi >n , g 2 ,n and g 3jn (with n > 1), only starts from the 0(rj) term. As anticipated, the fine 
structure of these charges, that is their dependence on the integer n, resides in the subleading 
Oijf) terms. The semiclassical results (14.291) provide a very accurate description of the exact 
spectrum of the quantized charge $3, see Fig. [31 For example, for N = 20 we find the following 
accuracy in each order of the ^-expansion of q^x, Eq. ( 14.29)) . 

77: +4.77%, rj 2 : -0.09%, r/ 3 : -0.0003%, (4.31) 

confirming the fast converge of the semiclassical expansion. 



4.4 Semiclassical anomalous dimensions 

Having determined the quantized values of the conserved charges, we can immediately find per- 
turbative semiclassical expansion of the roots of the transfer matrix. When substituted into the 
general expressions for the anomalous dimensions (13.36)) - (13.381) . it yields the trajectories in the 
vicinity of the upper boundary of the spectrum enumerated by the integer n. Note that at any 
loop level the argument of psi-functions contains only large leading order roots 5^ = O(N) so 
that one can find approximate formulas by expanding these -^-functions in the limit 77 <C 1, or 
rfe o/^l ^ 1- This consideration results in the following series in powers of rj, 



7o 



-3 In (V /2 r/ e" 7E ^ - 3r/(n + 1) - + 8n + 5n 2 ^j 



-^f^ + ^ + fn+fl+O^), (4.32) 
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N N N 

Figure 4: The exact eigenvalues of anomalous dimension matrix in the lowest three orders of 
perturbation theory for twist-three operators obtained from numerical solution of the Baxter 
equation and selected semiclassical trajectories for n = 0,1,4,7 at corresponding orders in 't 
Hooft coupling determined by Eqs. (14.321) and (I4.33p . respectively. 

for one loop, and analogously for two- and three-loop anomalous dimensions 

71 = -^ 2 7o - |C(3) + iv (v 7 ^ + 6 7 S 01 ) + \r, 2 (e(n + 1) ( 7 J 01 - 3) + vftr(3n + 2) + |) 

+ \rf (-9n + V3ir (f| + lOn + In 2 ) + 2 (f| + 8n + 5n 2 ) 7J 01 - 33n 2 - 60n - ±§i) + 0(n 4 ) , 

72 = Ao + f C(5) + |vr 2 C(3) + §77 (-2?C(3) - |v / 3vr 2 (vr + 2^)) (4.33) 
+ y ( - fvr 2 - 3( 7 i 0] ) 2 - v/Svr^ 01 + 18 7 l° ] + Sv^vr - 27C(3)(n + 1) 

+ |7r 2 [18(n + 1) - 6 7 J 0l (n + 1) - >/37r(3n + 2)]) + <D(rf) , 

where 7 i° ] = -3 In (3 1/2 r/e" 7E ). In Fig. 1 we compare the exact spectra at one, two and three 
loops with above semiclassical trajectories. In particular, for iV = 20, the accuracy of the two- 
loop asymptotic formula for 7l order-by-order in 77 is 

n: +0.37%, rj 2 : +0.11%, r/ 3 : +0.08%. (4.34) 

We observe that, in each order of perturbation theory, the anomalous dimensions given by 
(I4.32p and (14.331) receive corrections ~ ln(l/ry) growing logarithmically with the conformal spin 
r]~ l ~ N. Retaining the dominant contribution in each order of the ^-expansion we obtain 

1n,l=3 = -3 [g 2 \n7] + g A (-^TC 2 \nr] + §77 In 77) + g 6 (^ir 4 \nn -\n\nn + f^lnr/) 2 )] + . . . 

(4-35) 

where ellipses denote corrections suppressed by powers of n ~ 1/N. It is straightforward to verify 
that the same relation can be rewritten as 

1n,l=3 = -3T cusp (g) [\nn + lT cusp (g)r]\nn + §(r cusp (si)?7 \nn) 2 ] + . . . 

with r cusp (g) being the cusp anomalous dimension to three loops, Eq. (II. 8p . or equivalently 

7tv,l=3 = 3r cusp (s) In (N + §r cusp (s) In N) + . . . (4.36) 
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in a perfect agreement with (11.131) . 

The relation (14.361) suggests that the appearance of In 77 terms in the expansion of the anoma- 
lous dimension is an artifact of improper choice of the expansion parameter. The expansion can 
be recast in a form free from the large terms by recalling that in higher orders of perturbation 
theory the total conformal spin gets shifted by the anomalous dimension of the Wilson operator 
with the corresponding scaling dimension [49] 



Jo 



J 



N + U + h N , L {9) ■ 



This implies that in higher orders the anomalous dimension 7jv,L is actually a function of the 
renormalized spin j [28J, i.e., 

7tv,l = h(N + \ lN>L ) . (4.37) 
This relation can be easily inverted and the first few orders of its perturbative series read 
f(N) = j(N) - i(7 2 (A0)' + ^(7 3 (A0)" + <?(7 4 (A0) 



9 7o(?7) + 9 li{v) + \v lo(vho'(v) 



(4.38) 



+9" 



72M + W (70 (vhMY + W^o(v)) 2 io'(v) 

+ \V 4 (|(7o(r/))V(r/) + 7o(r/)(7o'(r/)) 2 ) 1 + 0(g s 



with the expansion coefficient expressed in terms of 7„ given in Eqs. (14.321) - (I4.33p . 

Since the auxiliary charge q<i depends on the renormalized conformal Casimir, we may re- 
express anomalous dimensions in terms of the parameter v = l/(— ^(fiO) 1 ^ 2 related to it rather 
than rj. A straightforward rearrangement of the semiclassical series yields 



where 



7o 



7i 



72 



7tv,l=3 = 9 2 lo(v) + g%(v) + g%(v) + 0(g 8 ) , 

-3 ln(3 1/2 w e" 7E ) - 3 (n + §) v - \ (5n 2 + 5n + f ) v 2 

-| (n + |) (29n 2 + 29n + ±f ) v 3 + 0(v A ) , 

-T2* 2 lo - |C(3) + + § (v^tt (n + |) + |) v 2 

+ \ (v^vr (7n 2 + 7n + g)-18 (n + l))v 3 + 0(v*) 

^ Ao + f C(5) + ^ 2 C(3) - i ^3- 



7T V 



(4.39) 

(4.40) 
(4.41) 

(4.42) 



-f (fvr 2 + v^vr 3 ( n+ ±))v 2 + 0(v 3 



As it was mentioned above, is a function of the bare total spin jo such that the N- 

dependence of /z,(iV) originates solely from the one of the bare charge q2^(N). Therefore, re- 
expanding fL(N) in the v = l/(— g 2 ,o) 1 ^ 2 power series we finally obtain 



h=s(N) = g 2 Uv ) + g 4 fi{v ) + g 6 Mv ) + 0(g* 



(4.43) 
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with v = [(N + §)(iV + |) + W 1 ' 2 and 



/ = -31n(3 1 / 2 Wo e-^) - 3 (n + §) « - ± (5(n + ±) 2 + g) u 2 (4.44) 

-|(n+|) (29(n + i) 2 + f) ^ 3 + OK 4 ), 

/i = -^vr 2 /o - JC(3) + iVSnvo + f (V^n + §) + §) «g (4-45) 

+ \ (v/Stt (7(n + |) 2 + f ) - 9(n + |)) ^ 3 + 0(^ 4 ) , 
/a = #^ 4 /o + f C(5) + |vr 2 C(3) - frfo?v Q (4.46) 

We recall that nonnegative integer n enumerates the trajectories close to the upper boundary of 
the band (jl.6p with the latter corresponding to n = 0. We observe that, in agreement with our 
expectations, the expansion coefficients in front of power suppressed corrections on the right-hand 
side of (14.441) - 1 14. 46ft are free from logarithmically enhanced contributions. 

5 Discussion and conclusion 

In this paper we analyzed the fine structure of the spectrum of the anomalous dimensions of 
Wilson operators in the autonomous SL(2) sector of the maximally supersymmetric gauge theory 
in the few lowest orders of the perturbative series in 't Hooft coupling. We identified two sets of 
trajectories corresponding to excited states in the vicinity of the upper and lower boundary of 
the band 111.61) which define two distinct analytical continuations for the anomalous dimensions. 
Both of them are parameterized by sets of integers which are related to each other through a 
set of linear relations. We focused in great detail on twist-three operators which, while being 
simple enough for analytical treatment, contain all salient features involved in more general 
considerations of twist-L spectra. 

The integrability of the dilatation operator played a key role in our study. We have applied the 
method of the Baxter equation, generalized to all orders of perturbation theory, for the eigenvalues 
of the Baxter polynomial which encodes the spectrum of anomalous dimensions of twist operators. 
Though the equation is only asymptotic in its nature, i.e., it does not implement wrapping 
effects of the putative long-range spin chain underlying the dilatation operator of the M = 4 
super Yang-Mills theory, this difficulty did not show up in the present three-loop considerations. 
In the past, the Baxter equation was efficiently used to study asymptotic regimes in anomalous 
dimensions at one loop order and was presently extended to multiloop analyses. The main focus 
of our consideration was the construction of a systematic expansion in inverse powers of the 
Lorentz spin which define corrections to the logarithmic Sudakov scaling. While the latter is 
universal to all gauge theories, the preasymptotic effects are not and strongly depend on whether 
the underlying dynamical system is integrable or not. We developed efficient procedures for 
systematic studies of anomalous dimensions in higher orders of perturbation theory generalizing 
techniques relying on the Baxter equation used in the past to study short-range noncompact spin 
chains in the vicinity of the lower and upper band of the spectrum. Due to different scaling of 
the conserved charges parameterizing the all-order auxiliary transfer matrix entering the Baxter 
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equation depending on the location within the energy band, we relied on either asymptotic 
or semiclassical techniques. Our analysis yielded a set of analytical formulas which accurately 
describe the spectra of anomalous dimension exhibiting very fast convergence of asymptotic and 
semiclassical series within the domains of their validity. 

The trajectories, being eigenvalues of certain hermitian long-range Hamiltonian of the puta- 
tive AdS/CFT long-range spin chain with different quantum numbes, do not intersect each other. 
However, their distribution density does change as one varies the value of the 't Hooft coupling 
constant and other global parameters of Wilson operators. This raises the question of precise 
identification of stringy configurations dual to the excited trajectories at strong coupling and 
corresponding dual observables. It is well known that the lowest trajectory corresponds to the 
double folded closed string extended to the boundary of the anti-de Sitter space and spinning with 
large angular momentum. The excited trajectories correspond to the spiky string configuration 
with its spikes in the vicinity of the boundary. Along these lines one would gain understanding 
of the density of states inside the band for large spins at strong coupling. 

Another problem which emerges from our consideration is the study of anomalous dimensions 
of twist- L operators as one moves into the interior of the band starting from the low boundary. 
Along the way, one hits a special set of trajectories (11. lip which scale at large Lorentz spin as 
(K+2) InN. From the point of view of the spectral curve underlying the magnet this corresponds 
to opening double points on the Riemann surface of the auxiliary problem as one goes upwards the 
band of eigenvalues. This continuous deformation of the density of Bethe roots can be described 
through the Whitham flow of energy eigenvalues with respect to the change of the integrals of 
motion. These issues will be analyzed in forthcoming publications. 
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